Zipf's law in Nuclear Multifragmentation and Percolation Theory 
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We investigate tlie average sizes of the n largest fragments in nuclear multifragmentation events 
near the critical point of the nuclear matter phase diagram. We perform analytic calculations 
employing Poisson statistics as well as Monte Carlo simulations of the percolation type. We find 
that previous claims of manifestations of Zipf's Law in the rank-ordered fragment size distributions 
are not born out in our result, neither in finite nor infinite systems. Instead, we find that Zipf- 
Mandelbrot distributions are needed to describe the results, and we show how one can derive them 
in the infinite size limit. However, we agree with previous authors that the investigation of rank- 
ordered fragment size distributions is an alternative way to look for the critical point in the nuclear 
matter diagram. 
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I. INTRODUCTION 

A central goal of nuclear physics is to study the phase 
diagram of nuclear and quark-gluon matter. Of special 
interest are two phase transitions, both of which are ex- 
pected to be first-order transitions terminating at a crit- 
ical point in the — T plane. At high temperatures 
(T w 200 MeV) and/or densities QCD tells us that nu- 
cleons are no longer bound, but dissolve into assymp- 
totically free quarks and gluons possibly with an 
accompanying restoration of chiral symmetry The 
SPS-CERN and RHIC-BNL experiments have found ev- 
idence of this phase transition, and that is might be of 
first-order (or possibly a cross-over) [3|, |4| . Theoretically, 
lattice-QCD calculations suggest that the nuclear matter 
phase diagram contains a critical point at high tempera- 
ture and density ^, d, Q • One of the central goals of the 
future FAIR facility at the German GSI is to study this 
critical point. 

The second transition is under far better control, both 
experimentally and theoretically, and involves the tran- 
sition of the low-temperature nuclear Fermi liquid into 
a gas of nucleons and small fragments. Experimentally 
[i,H[i3,[ni,[Il,IIl,[Il, one finds a critical point at which 
this first-order transition terminates. This critical point 
is located at temperature and density [l^ [l^ [iTj 



Tc = 8.3 ± 0.2 MeV, p^/po = 0.35 ± 0.1 (1) 

In addition, through comparisons with phase transi- 
tion theories, the set of critical exponents determining 
the universality class of the phase transition has been 
determined [l5| . The experimental values of the expo- 
nents, extracted in a somewhat model dependent way, 
are in good agreement with the universality class of per- 
colation [lllia,!!^,!!!!!!!!!, provided that appropriate 
finite size corrections are included [Ts, 24, 25]. In passing 
we note that similar fragmentation patterns pointing to 
the existence of phase transitions can also be found in the 



disintegration of molecules, such as buckyballs [26|, |27|. 
With the advent of the very high-powered free-electron 
lasers at DESY in Hamburg and at SLAG in Stanford a 
fresh look at these molecular fragmentation experiments 
will be possible, and first experiments have been pro- 
posed [28(. 

The great majority of analyses of critical phenomena 
in nuclear fragmentation have considered the shape of the 
mass distribution, which behaves as a power law in the 
critical region. In the present article we wish to address 
whether studying the behavior of the ranked fragment 
sizes, i.e., the sizes of the largest fragment, second largest, 
and so on, also contains information regarding critical 
behavior. As a second goal, we also investigate the ap- 
plicability Zipf's law, an empirical expression sometimes 
linked to self-ordered criticality which has recently been 
applied to nuclear multifragmentation. 



II. ZIPF'S LAW IN NUCLEAR 
MULTIFRAGMENTION 

When Watanabe ^ and Ma et al. [H [IH studied 
the average size of the n^^ largest fragment (An) as a 
function of n they found that the data is described by 
Zipf's law [32] which was first found in linguistics and 
is named after the linguist George Zipf). Zipf's law is 
empirical, and states that the most frequent word in a 
given language appears twice as often as the second most 
frequent word, three times as often as the third most 
frequent and so on, i.e. 



F{r) oc r 



(2) 



where F{r) is the frequency of the r-th most frequent 
word. A more general form. 



-A 



F{r) cx r 

is also often refered to as Zipf's Law. 



(3) 
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Since it was first formulated, many examples outside 
of linguistics have shown the same behavior. It should 
also be noted that the Pareto distribution [s^ (named 
after the economist Vilfredo Pareto) contains Zipf's law 
in a different formulation and was also found in a broad 
range of statistical data. 

Ma et al. claim that since Zipf's law has been observed 
in very different fields this is a reflection of self-ordered 
criticality. Therefore, in [s^l Zipf's law in its extended 
form ^ was fitted to experimental data for the aver- 
age size of the n^^ largest fragment, see Figs. 23 and 24 
in [30| . They found that for increasing excitation ener- 
gies A decreases and passes through A = 1 at some point, 
i.e. Zipf's law in its original version It is hoped that 
this passing through A = 1 is an indication for the critical 
point. 

III. n-TH LARGEST CLUSTERS AND 
AVERAGE SIZES 

If we assume that the production of any two clusters is 
largely independent, then the probability that a cluster 
of size A is the n*''-largest cluster is given by 

P„(A) = ^P>„_,(^).k(>^) , (4) 

i=0 

where 

m — 1 

p>„^{A) = 1 - PkiA) (5) 

fc=0 

is the probability to have at least m clusters of size A. 

Here, Pk{A) is the probability to have k clusters of size 
A and Pi{> A) is the probability that there are i clusters 
of size larger than A. For example, the probability that a 
cluster is the largest in a given fragmentation event is the 
product of the probability p>i{A) that there is at least 
one cluster of size A present and the probability Pq{> A) 
that there are zero clusters of size larger than A 

Pi{A) = p>i{A) ■ poi> A) 

= il-PoiA))-poi>A) (6) 

The probability that a cluster is the second largest is 
given by a sum of two terms. The first term is the prob- 
ability p>2 (A) that there are two or more clusters of size 
A present multiplied with the probability po{> A) that 
there are no clusters of size larger than A. The second 
term is the probability p>i{A) that there is at least one 
cluster of size A present multiplied with the probability 
Pi{> A) that there is exactly one cluster of size larger 
than A present. This yields 

P2{A) = p>2{A) ■ po{> A) + p>i{A) ■ pi{> A) 
= {1 - po{A) - p,{A)) ■ po{> A) 

+{1 - Po{A)) ■ p,{> A) (7) 



for the probability that a cluster is the second largest in 
a given fragmentation event. From the probability P„th 
that a cluster of size A is the n*^ largest we can calculate 
the average size of the n-th largest cluster 

V 

(A„) = ^A.P„(A) . (8) 

Here, Pn{A) depends on the probability distribution pi 
that is chosen for the underlying physical system. Since 
we assume that the production of any two clusters is 
largely independent from one another, the probability 
distribution can be approximated by the Poisson distri- 
bution 

p,(A) = l(7V(A)f e-^(^) , (9) 

where N(A) is the appropriate cluster size distribution, 
i.e. the average of clusters of size A, and may depend on 
additional parameters like the excitation energy of the 
reaction. 

For a system close to the critical point, the cluster size 
distribution follows a scaling function of the form 

n(A,e) = aA-VM") , (10) 

where e is the fractional deviation of the control parame- 
ter (for example the temperature T, i.e. e = [T — Tc)/Tc) 
from the critical value. The scaling function / is equal 
to 1 at the critical point (/(O) = 1). The critical ex- 
ponents T and a determine the universality class of the 
phase transition, and the normalization constant a is de- 
fined by the condition that all nucleons belong to some 
cluster 

V 

Y^A-n{A,e)^V , (11) 

A=l 

where V is the total number of nucleons in an event. 

In [s^l and [35| the behavior of the average sizes of the 
n-th largest cluster at the critical point (e = 0) was stud- 
ied assuming that the cluster sizes can be described by a 
Poissonian probablity distribution. Here, we summarize 
those results and consider the critical point where the 
cluster size distribution is a pure power law 

N{A)^n{A,{))=aA-^ , (12) 

which with equation ^ can be used to calculate the av- 
erage size of the n-th largest cluster according to equa- 
tions Q and ([HI) and we will consider 2 < r < 3 in the 
following. The normalization constant is given by 

V 

a = V/Y, ^-(1-^) = V/Hv,i-r , (13) 

A=l 

where Hn^m = the n*^ harmonic number of 

order m. For the probability to have i-clusters larger 
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than A, 

V V 

N{>A) = J2 N{k)=a J2 (14) 

k=A+l k=A+l 

= a [C(t, 1 + A)- C(t, 1 + V)] 

has to be used with equation where 
(^(s, g) = ^j.^g(A: + g)^'' is the generaUzed Riemann 
function. 

For large systems, one can also replace the sum in equa- 
tions to (fTSl) by an integral. The constant a is then 
defined by 



dAaA-^^-^^ = V 



(15) 



and is given by 

a = V{r-2)[A;J-'^^V-(^-'Y ' 

For the average number of fragments larger then A we 
obtain 



7j{A) = [ dAaA-^ 

J A 

a 



T - 1 



The probability that a cluster is the n^^ largest is 

„n-l 



Pn{A) = 



77" 



(n-1)! 



(16) 



(17) 



(18) 



and therefore the average size of the n largest cluster 
is 

rv n-l 

(A„)-/ dAAaA-^^^—e-"^. (19) 
JA„.„ (n-l)! 

We can rewrite this integral by integrating over ri{A) with 

V ^ 

and get 



Ai,)^(^-^ + V-i^-A'^-' , (20) 



r- 1 



drj [ T] + a 

V ^-1 

„n-l 



X e 



(n-l) 



(21) 



where rj{V) = and rj-^sx = ??(^min)- For average size of 
the n + 1-th largest fragment in an infinite system this 
leads to 



T- 1 



drj e 



n"-—- 
-n 1 

(n-l)! 



(22) 



which we can express in terms of the gamma function as 

^r(n-^) 



= 



T - 1 



(n-l)! 

^r(n-^ 
r(n) 



(23) 



where in the last step we made use of the fact that r(n) — 
(n-l)!. 

Hence, in the case of an infinite system for r = 2 the 
average size of the n*'' largest cluster is 



r(n 



r(n) 



^ = a(n-l)-i , 



(24) 



where we used the identity r(a; -f 1) = a;r(a;). Therefore 
for T = 2 the average size of the n"^ largest cluster in an 
infinite system is not described by Zipf 's law, but by the 
Zipf-Mandelbrot distribution [H, [s^l 



F{r) ^ c{r + ky 



(25) 



with k — —1 and A = 1. Fig. [T] shows {An) calculated 
with equation ([5]) , divided by the normalization constant 
a for different system sizes V . For n = 1 and n = 2 the 
finite size effects are very noticable, however for n > 3 the 
{An) match the Zipf-Mandelbrot distribution. It should 
also be noted that Zipf 's law with A = 1 does not describe 
{An) very well. 



2 

10 

5 



"a 

< 

V 





A V=10^ ■ 




• v=io'* 




■ V=10^ - 






1 \^ 


Zipf 















FIG. 1: Symbols indicate {An) for t = 2 and different sys- 
tem sizes calculated with equation , the solid line indicates 
{A^) from equation (|24p . the dashed line indicates Zipf's 
law ©. 

Using the ansatz described above we determine {An) 
for different values 2 < r < 3 and test how well they are 
described by Zipf's law and the Zipf-Mandelbrot distri- 
bution. 

Fig. [2] shows the extracted fits of the parameter A for 
Zipf's law for different values r and different system sizes 
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FIG. 2: Top; A as a function of r as fitted to Zipf's law Q 
for different system sizes. Bottom: Error of the fit 




V . The value of A for the finite systems shows a decrease 
from A(t = 2) = 1.25 to A(r = 3) w 0.75, whereas for 
the Hmit of an infinitely large system there is a stronger 
dependence on r. The errors of the fit are shown in the 
lower part of the figure, demonstrating that the fits for 
the finite systems are a good approximation, while the 
error for the infinite system is substantially higher. For 
the finite systems, A is relatively close to unity for the 
small system with V — Kfi around r — 2.5. However, for 
example A(t = 3) ~ 0.75 which would indicate that the 
statement of [13] that A = 1 as found in the experimental 
data is more coincidental than profound. Fig. [2] shows 
that the value of A at the critical point depends on r and 
the system size V and that A = 1 at the critical point 
only holds true for certain sets of parameters. 

Fig. [3] shows the extracted fits of the parameter A 
for the Zipf-Mandelbrot distribution for different r and 
different system sizes V . The value of A for the fi- 
nite systems shows a decrease from A(r = 2) w 1.2 to 
A(t = 3) w 0.5, while for the larger system the decrease 
is not as strong. The difference between the finite sys- 
tems and the infinite system is not as strong as for the 
fit to Zipf's law. The errors of the fit are shown in the 
lower part of the figure showing that the fits for the finite 
systems are indeed a much better approximation of the 
{An) than Zipf's law, and the fit for the infinite system 
yields the best approximation overall. The parameter k 
shown in Fig. [3] (middle) shows some finite size depen- 
dency for small r, however the system size dependence 
decreases with increasing t. 

From fits to the Zipf-Mandelbrot distribution for 2 < 
r < 3 one finds a systematic dependence of A on the 



FIG. 3: A as a function of r as fitted to the Zipf-Mandelbrot 
distribution (|25|l for different system sizes. Middle: Extracted 
parameter k. Bottom: Error of the fit 



critical exponent r of the form 

1 



A 



r - 1 



(26) 



which holds true almost exaclty for the infinite system 
and is a good approximation for small systems. 

To summarize the findings of this section: Zipf's law 
should not be applied to nulcear fragmentation, but 
rather the more general form of the Zipf-Mandelbrot dis- 
tribution. And A = 1 as found in the experimental data 
is more coincidental than a measure of a critical exponent 
or an indicator of the self-similar behavior. 



IV. PERCOLATION THEORY AND ZIPF'S LAW 

The ansatz described in IIIII can only be applied for 
systems at (or close to) the critical point since informa- 
tion about the scaling function / is needed for e 7^ 0. It 
has been discussed in [l^l that bond percolation can be 
used to describe multi-fragmentation reactions. In the 
following we will use three dimensional bond percolation 
model on a simple cubic lattice 0, [111 . The probability 
to form a bond between two neighboring lattice sites is 

Pbond- 

The critical exponents of the percolation model depend 
on the system size and the dimensionality of the lattice. 
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For an infinitely large three dimensional lattice r = 2.18 
and a = 0.45. The critical bond probability does not 
only depend on the dimensionality of a lattice, but also 
on its particular topology. For an infinitely large simple 
cubic lattice the critical probability is Pc = 0.249. 

In the following, we will consider two different system 
sizes: a small lattice with 6x6x6 = 216 sites (corre- 
sponding to system sizes found in nuclear fi'agmentation 
reactions) and a large lattice with 100 x 100 x 100 — 10^ 
sites (reducing the finite size effects considerably). 




1 2 5 10 1 2 5 10 




1 2 5 10 1 2 5 10 




FIG. 4: Average sizes of the n-th largest cluster as a function 
of n for the lattice with 216 sites for different bond probabil- 
ities pbond. The dashed lines are fits for Zipf's law Q, the 
solid lines are fits for the Zipf-Mandelbrot distribution (|25[). 

Fig.[4]shows the average sizes of the n-th largest cluster 
as a function of n for the lattice with 216 sites for dif- 
ferent bond probabilities pbond- For this system size we 
effectively have r w 2.14 at the critical point « 0.253. 
The value for pc was chosen by determining the value 
of Pbond where the cluster size distributions N{A) were 
fitted best by the power law distribution (fT2)) . 

The results for {An) from percolation theory are simi- 
lar to the experimental data shown in Fig. 23 from [30| . 
The small Pbond correspond to high excitation energies 



and vice versa. For our results from percolation we see 
that the fits for Zipf's law (and Zipf-Mandelbrot) work 
very well below Pc (that corresponds to excitation en- 
ergies above the critical point). But as one crosses the 
critical point to pbond > Pc (corresponding to excitation 
energies above the critical point) we observe that Zipf's 
law (and eventually Zipf-Mandelbrot) does not work if 
(Ai) is included. 

We fit Zipf's law ([3]) to (An) from percolation and de- 
termined the parameter A. Since we use a least squares 
fitting procedure we use two different methods to obtain 
a fit. The first method (denoted "Fit 1") is equivalent 
to what is done for Fig. 24 in i.e. fitting the data 
to equation Q with a least square fit. However, this 
means that for a power law fit the first two data points 
will dominate the result for the fit. As a second method 
(denoted "Fit 2" ) we first take the logarithm of both axis 
and then fit the (An) to a straight line to determine A, 
hence all the (An) contribute more or less equally to the 
fit result. 

We already mentioned that Zipf's law (Zipf- 
Mandelbrot) does not work if the average size of the 
largest fragement (Ai) is included above Pc- Therefore 
we perform the two fit methods once including all the 
fragments and once excluding {Ai) from the fit. This 
way we have four different ways to determine a value for 
A and the results are shown in Fig. [5l for both the small 
and the large lattice. 

For the small system we see a similar behavior as for 
the experimental data from [s^l and A(pbond) crosses A — 
1 at or very close to Pc for both system sizes investigated. 
The fits are all very similar below Pc as the first fragment 
hardly deviates from Zipf's law. 

In case of the large system (10^ sites), and for the fit 
including (Ai), a quick rise of A can be observed with 
increasing pbond due to the domination of (Ai). For the 
fit excluding (Ai) we see a distinct peak of A(pbond) close 
to the critical point. 

Although some maximum is visible for the case of the 
small system (216 sites), the maximum is very broad and 
not particular near the critical point and can therefore 
not be taken as an indication of the critical point for 
these small systems. 

If we fit (An) with the Zipf-Mandelbrot distribu- 
tion ([25|) . we get a similar behavior for A. Further, for 
the case where we exclude the largest fragment from the 
fit, the parameter k also indicates the critical point in 
the case of the large lattice (10^ sites). Here we observe 
that fc « for pbond < Pc and a sharp drop to fc « —1 for 
Pbond > Pc, i-e. the dependence of k can be approximated 
by 

fc(Pbond) ~ -S(Pbond -Pc) ■ (27) 

This means that we recover Zipf's law if we neglect the 
largest fragment above p^. However, when turning to the 
small lattice (216 sites), the sharp drop of k at pc turns 
into an extremely smooth decrease over a range in the 
bond probability of Apbond ~ 0.2. 
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FIG. 5: Fitted values for A using the different fits described 
in the text. Solid symbols are for fits including the first frag- 
ment, open symbols are for fits omitting (An). Upper figure is 
for the small lattice (216 sites), lower figure is for the large lat- 
tice (10^ sites). Approximate values for pc are also indicated 
by a dashed line for both cases. 



Therefore, the crossing through A = 1 in the experi- 
mental data could indeed indicate the crossing through 
the critical point. However, we should keep in mind the 
results from sect. IIIII which clarified that the value of A 
at the critical point depends on r and V. Therefore, the 
fact that A = 1 close to the critical point cannot be taken 
as a unambigous signature for the critical point. 



V. PERCOLATION THEORY AND {A„) 

We have argued in the previous section that the Zipf 's 
law parameter A is not a rigorous signature for the critical 
point. In this section we want to investigate what can 
be learned from considering the average size of the n-th 
largest clusters (An). 



We turn to look at the ratio of the n-th largest average 
cluster sizes If there is a critical behavior 

that indicates the critical point we hope to find it by 
looking at the ratio of the (A„) for different n. 



V = 100 




0.25 

Pbond 

FIG. 6: Ratio of the average size of the n + 1-th largest frag- 
ment to average size of the n-th largest fragment as a function 
of Pbond for a large system with 10® sites. The result for n = 1 
is omitted. 

Fig. [6] shows the ratio (A„+i)/(A„) (with n — 2 . . .7) 
as a function of pbond for the large system (10^ sites). 
Close to the critical point all ratios show a significant 
minimum at the critical bond probability pc- The aver- 
age size of the largest fragment (Ai) has a very different 
behavior compared to (An) for n > 1 and is therefore 
omitted. 

Now we want to compare the results shown in Fig.[S]to 
the ansatz from section Hill for an infinite system. In sec- 
tion |TTT] two assumptions were made: at the critical point 
the cluster size distributions are discribed by a power 
law and the production of the clusters is uncorrelated 
and therefore Poissonian. 

For percolation theory we also expect to have a power 
law for the cluster size distributions at the critical point. 
Therefore, the only difference we would expect at the 
critical point compared to the ansatz from section IIIII 
for an infinite system besides finite size effects, would 
be due to correlations. These correletions would stem 
from the fact that the production of large clusters is not 
independent. For example, if a single very large cluster 
is produced whose size is of the order of the total sytem 
size no further large clusters can be produced in the same 
event. However, the production of the smaller clusters 
would not be affected. Hence, at the critical point we 
wouldn't expect the Poissonian ansatz to work for the 
ratio {An+i)/{An) for n = 1, but possibly for n larger 
than 2 or 3. 

Before comparing the percolation results to the ansatz 
from section IIIII for an infinite system we first want to 
make sure that finite size effects are negligible. Fig. [7] 
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shows the ratio {An+i) / {An) for different system sizes. It 




1 2 3 456789 



n 

FIG. 7: The ratio {An+i) / {An) for different finite system sizes 
V (symbols) using eqn. ((8)1 and for the infinte system (line) 
using eqn. ^ for r = 2.18. 

shows that finite size effects between the infinite system 
and the finite system sizes are negligible for n > 1 for 
the system sizes larger than V = 10^, and for n > 2 
are negligible for all shown system sizes. Therefore, the 
ansatz from section Hill for an infinite system can be used 
for the comparison to the percolation model if the ratio 
for rt = 1 (and rt = 2 for the smaller system) is excluded. 
However, it should be noted that for r = 2 the finite size 
effects are again substantial for n = 1 and non-negligible 
for n — 2 and n = 3, the finite size effects get negligible 
for n > 3, as is shown in Fig. [5] 




12 5 10 



n 

FIG. 8: The ratio {A„+i) / (An) for different finite system sizes 
V (symbols) using eqn. ([SJ and for the infinte system (line) 
using eqn. (|24p for t — 2. 

To compare the percolation result to the results from 
ansatz from section Hill that resulted in equation (P51) . 
Using the identity r(a; + 1) = a;r(a;) we can write the 



ratio as 




0.8 



0.1 0.15 0.2 0.25 0.3 0.35 0.4 

Pbond 

FIG. 9: Ratio of the average size of the n + 1-th largest frag- 
ment to average size of the n-th largest fragment as a function 
of Pbond normalized by the factor in equation (|28|l for a large 
system (10^ sites). The result for n = 1 is omitted. 

The ratio of the percolation result and the ansatz from 
section IIIII as a function of pbond is shown in Fig. [9l The 
minima all come down to 1 for n > 3 and one can con- 
clude that for the percolation model (at the critical point 
where the fragments are distributed according to a power 
law) the production of the n-th largest cluster with n > 2 
is not correlated for a large system. 



, Peril ~ 0.253 

V = 6 1 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 

Pbond 

FIG. 10: Ratio of the average size of the n + 1-th largest 
fragment to average size of the n-th largest fragment as a 
function of Pbond normalized by the factor in equation (|28|) for 
a small system (216 sites). The result for n = 1 is omitted. 
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When turning to the small system (216 sites) shown in 
Fig. [TUlhowever. the strong indication of the critical point 
by the location of the pronounced minima does not hold 
true anymore. We observe a broadening of the minima 
and a shift of the minima away from the critical point 
towards higher pbond- For n > 5 the minima still meet 
at some Pbond, but it is not located close to Pc ~ 0.253. 
This means that the percolation model and ansatz from 
section [ITT] agree for values of pbond where the cluster 
distribution is no longer a pure power law. 

VI. SUMMARY 

In contrast to previous claims, we have shown that a 
strict application of Zipf 's law to rank-ordered fragment 
size distributions near the critical point of the nuclear 
fragmentation phase transition is not warranted. This 
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tions. In order to obtain information on finite-size effects, 
we have performed these calculations over a wide variety 
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